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Abstract 

A central tenet of the new theory of gravity proposed by H. Yil- 
maz is the inclusion of a gravitational stress-energy tensor — t^ v along 
with the matter stress-energy tensor on the right hand side of 
the Einstein field equations. This change does not effect the New- 
tonian limit of the field equations since these terms are quadratic in 
potential gradients. From the Bianchi identities, however, important 
changes appear in any equations of motion consistent with these field 
equations. For matter described as a perfect fluid, and with Yilmaz's 
choice of signs when introducing these quadratic terms, we find that 
the Euler hydrodynamic equation in the Newtonian limit is modified 
to remove all gravitational forces. This allows, e.g., a solar system in 
which the Sun and the planets are permanently at rest, but does not 
explain how fluid bodies such as the Sun or Jupiter could form or be 
prevented from dispersing. 

PACS 04.20.Cv - Fundamental problems and general formalism. 

PACS 04.50.+h - . . . other theories of gravitation. 

PACS 04.25.— g - Approximation methods; equations of motion. 



1 Introduction 

Recent publications by Yilmaz ffl, |2[ claim that General Relativity (GR) is 
inconsistent and that its Newtonian limit is unsatisfactory. Because these 
claims have resulted in some widely distributed notice [Q, ||, and because 
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they may be having an invalid influence on experimental proposals ||, it 
seems appropriate to restate that they are incorrect. But the straightforward 
and long established arguments (or more rigorous recent work JjJ) showing 
GR to have a good Newtonian limit leave little room for including Newtonian 
stresses in still another way in the field equations, as Yilmaz has proposed. 
We will see here that the Yilmaz theory makes changes to GR that are so 
profound they remove all gravitational interactions at the Newtonian level 
in any exact solution of the field equations of the modified theory. 

In notation that differs in detail from that used by Yilmaz in his papers, 
but adopts the MTW conventions S, the Einstein field equations are 

where all components of T^ v have units of energy density (i.e., pressure) 
and this tensor describes the stress-energy of all non-gravitational fields. In 
the Yilmaz theory this equation is modified to read 

GT = (8vrG/c 4 )(T^ - IT) , (1.2) 

where in every instance t^ u as used here differs in sign from the quantity 
Yilmaz represents by this symbol. The exact definition of the Yilmaz grav- 
itational stress-energy tensor — t^ v is in some cases unclear, but explicit 
expressions have been given in situations that include the Newtonian limit. 
In suitable (harmonic) coordinate systems chosen by Yilmaz the definition 
becomes = Un where is the Einstein stress-energy pseudotensor; 
my meaning of t^ u avoids a minus sign in this definition, as well as in the 
relationship ( |4.3| ) below of t^ u to the Newtonian gravitational stresses. 

2 Relativistic Matter 

To compare the GR and Yilmaz theories in the Newtonian limit one needs a 
theory of matter so there is something to exhibit gravitational interactions. 
I will take for a theory of matter the relativistic hydrodynamics of a perfect 
fluid. [Newtonian hydrodynamics with gravitational forces is summarized 
in appendix |A].] This is a suitable test case as it avoids dealing with the 
singularities of point particles, yet includes some nongravitational (pressure 
gradient) forces which are needed in the Newtonian theory for static solu- 
tions to exist. In this case one takes the stress-energy tensor of matter to 
be 

= pc 2 u»u v + P(g^ + u»u u ) . (2.1) 
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The fluid 4-velocity is normalized by the condition 

9^vru u = -l (2.2) 
so that in quasi-rectangular coordinates the 4-velocity can be written 

«" = (l,v/c) 7 (2.3) 

where v k = dx k /dt is the coordinate velocity of the fluid particles, x° = ct 
is the fourth coordinate in tensor expressions, and 7 will be determined by 
the normalization ( |2.2| ). All components and u 11 are dimensionless, as 7 
will then be. 

In the Yilmaz theory the hydrodynamic equations have gravitational in- 
fluences that differ from those in GR. In this case of a matter theory with 
only four independent fields (p and the three v k with P(p) given by an 
equation of state) the Bianchi identities which enforce local energy and mo- 
mentum conservation constrain the matter equations of motion sufficiently 
to determine them completely. The contracted Bianchi identities G^ v - V = 
from equation ( |1.2| ) imply {T^ u — t^ v )- u = in the Yilmaz theory. 

When split into components parallel and perpendicular to the fluid 4- 
velocity the local conservations equations {T^ v — t^ u )- v = are an energy 
equation 

{p?u* + v v r*)w = -Pu v . v + ti a 'uKv , (2-4) 
from —u^{T^ v — t^ v )^ v = and the generalized Euler equation 

(pc 2 + P)u», v u v = -{sT + u»u»)(P v - t u a . a ) (2.5) 

from {5£ + u^u u )(T l ' a — t l ' a )- a = 0. The terms involving t^ u are absent 
for the Einstein theory of gravity. In equation ( |2.5| } this allows geodesic 
motion u fl - u u l/ = in GR when the pressure gradient forces on the right 
are also absent. Although the t^ v terms, like the pressure terms, should be 
negligible in the Newtonian limit in the gravitational field equations (|1.2| ) 
and also negligible in the equation of continuity ( p.4[ ) they may, like the 
pressure gradient terms, be important in the Newtonian limit of the Euler 
equation ( |2.5| ). We will find that they cancel the terms from pc^u^-yU 11 that 
in Einstein's theory give the Newtonian gravitational force term in the Euler 
equation. 

The generalized Euler equation ( |2.5D is not optional or modifiable in the 
Yilmaz theory. If a t^ v term is inserted in the Einstein equations, the Bianchi 
identities show that it must appear in the hydrodynamic equations of motion 
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satisfied in any exact solution. Thus to leading order it must be satisfied 
by any motion that approximates an exact solution of equations ( |1.2| ) with 
a hydrodynamic stress-energy tensor. This argument makes no use of 
any additional equations that may be part of the Yilmaz theory, as we have 
to this point made use only of the field equations (fj~^) and the assumption 
that the theory is applicable to selfgravitating systems of fluid objects with 



the stress-energy tensor (2.1 



3 An Exact Solution 

Yilmaz has given an exact solution of his field equations g ||] which I derive 
here to verify that, in spite of changes in notation and sign conventions, I am 
stating equations from his theory correctly. When the consequences fl2.5| ) 
of the field equations are considered in this example one sees a failure to 
reproduce Newtonian gravitational effects. The metric Yilmaz proposes is 

ds 2 = -e 2 *d(ct) 2 + e- 2 *(dx 2 + dy 2 + dz 2 ) . (3.1) 

For this choice of metric and coordinates Yilmaz evaluates his gravitational 
stress-energy tensor — t^ u to give 

V = \(c a /AkG)[<S>.^ - ±5fc a & a ] (3.2) 

with A = 1. We have inserted the tag A here so that by setting A = in later 
equations we can recover the corresponding equation in Einstein's theory. 

I now use the computer algebra aid GRTensor II || [l0| to construct the 
Yilmaz field equations in this case. The result is 

2$ kk + (A - 1)* fc * fc + (A + 3)e- 4 *(cf>, ) 2 
(87rG/c 2 )[p 7 2 + (P/c 2 )( e - 2 *-7 2 )] , (3.3) 

-2$ J0 - 2(A - 1)$^, 

(SttG/c^vHp + P/c 2 ) , (3.4) 

+<S y [2* 1 oo + (A-5)(* )0 ) 2 ] 
(8vrG/c 4 )[7 2 ^^( / o + P/c 2 ) + ^Pe 2 *] . (3.5) 

The left hand sides of these equations are the gravitational quantities Y^ v = 
G^ v + (87rG/c 4 )t /iiy while all the matter terms are collected on the right. 
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To find a static solution we set A = 1 (Yilmaz theory), drop all time 
derivatives, set ifl = 0, and from the normalization condition (2.2) get 7 2 = 
g-2* rp^g resu iti n g equations are 

(3.6) 
(3.7) 
(3.8) 

To construct solutions to these equations one chooses any desired density 
distribution p e s(x, y, z) representing one or more fluid bodies and then solves 
the Poisson equation U^k = ^Gp e s- Then the choice $> = U/c 2 reduces 
equation (|3.6|) to the form p e g = pe~ 2U ^° 2 which serves to display the in- 



yOO . 


2$,fcfc = 


{SnG/c 2 )pe- 2 * 


Y 0j : 


= 


o , 


yii . 


= 


{MG/c^Pe 2 * 



variant density p = u^u u T^ v /c 2 in terms of the given effective density. To 
solve equation ( |3.8| ) one must set P = 0. 

The existence of this exact solution is disastrous for the Yilmaz theory. 
Any equation of motion for fluid matter adjoined to the Yilmaz field equa- 
tions (|1.2| ) cannot differ from the generalized Euler equation ( |2.5| ) since that 
equation is just a combination of the field equations and their derivatives. 
Thus this static solution with P = satisfies also the fluid equations of 
motion (|2.5| ) in any consistent theory incorporating the Yilmaz field equa- 
tions (1.2). One such solution could represent two fluid planets (Jupiter 



and Saturn) placed near a fluid Sun with none of these bodies in motion. 
No gravitational forces act to accelerate any of these bodies, no mechanical 
(pressure) forces are needed to balance gravity, and the only gravitational 
effect is the renormalization between effective density p e ff and local invariant 
density p. 

The Einstein equations under these same static assumptions have no 
exact fluid solution in agreement with the Newtonian impossibility of a mo- 
tionless solar system. The usual Oppenheimer-Volkoff static fluid solutions 
(spherical symmetry only) are obtained by allowing —goo an d g xx to differ 
and do require pressure gradients to support the matter against its own 
gravitational fields. See, for example, ]8|, §23.7]. To see that allowing addi- 
tional independent components in g^ v does not help in the Yilmaz theory 
we treat the Newtonian limit below without assuming static configurations 
nor restrictions beyond harmonic coordinate conditions on the metric com- 
ponents. 
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4 Newtonian limit 



In the Newtonian limit one assumes that velocities are small, \v/c\ <C 1, 
that gravitational potentials g^ v are near their Minkowski values, and that 
pressures or other mechanical stresses are negligible compared to energy 
densities \P\ <C pc 2 . By arguments known since Einstein and Hilbert, one 
finds that only one combination of the g^ u is generated by the dominant 
term in T^ u which is pc 2 vPu°. (See, e.g., [||, §18.4].) Neglecting terms in 
g^ u smaller than this leading term by factors of v/c or <E> one finds that the 
metric is the same as equation (3.1) in which one may set e ±2 * = 1 db 2<£. 



This metric will be time dependent with (J?^ = d^/dct smaller than by 
a factor of v/c where v is the velocity of the matter sources. [Consider, for 
example, the Newtonian potential U = —GM/\r — vi| of a moving point 
mass.] The equality of —500 and g xx is deduced rather than assumed, and 
holds only as a first approximation. These arguments hold equally well for 
the Yilmaz theory or for the Einstein theory since t^ u is of order (d<&) 2 and 
thus smaller than the leading <9 2< 3? terms since |3>| <C 1. Thus both the 
Yilmaz theory and the Einstein theory lead to 

2$ fcJfe = (8vrG/c 2 )p (4.1) 

from the G 00 equation, while all terms in the other field equations are neg- 
ligible in comparison. The solution of this equation is 

$ = U/c 2 (4.2) 



where U is the Newtonian potential satisfying equation ( |A.l| ) for the same 
matter distribution p(x,y,z,t). 

Since both the Einstein and Yilmaz field equations give the same metric 
in the Newtonian limit, it had long been assumed that the Yilmaz theory 
had a satisfactory Newtonian limit. But Yilmaz M has recently suggested 
that the theories differ in this limit, prompting this study which localizes 
the difference to the equations of motion required by the Bianchi identities. 



Thus we proceed to study approximations to equations fl2.4|) and d2.5|) . 

The reduction of the relativistic continuity equation fl2.4j) is relatively 
straightforward using the metric (|3.1j) . The terms on the right by which 
restmass is changed by mechanical or gravitational work are unimportant 
for Newtonian strength forces, and in the mass current vector pc 2 u^ + u u t L ' IJ ' 
the gravitational binding energy contribution to the rest mass from is also 
small. This leaves (pcu^).n = in which the approximation = (l,w J /c) 
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will suffice. All terms involving <9<3? are smaller by a factor $ than those in 
the Newtonian equation and thus negligible, so equation ( |A.2| ) results. 

The relativistic Euler equation ( |2.5| ) similarly simplifies. (The p = 
equation is negligible.) The left hand side with u° = 1 and u l = v i /c 
reads pc 2 u l ^ u u u = p[(dv l /dt) + v k v\ k + T l ooc 2 ] and the right hand side is 
—Pi + t tk jk . But when one makes the identification <I> = U/c 2 from the 
field equations the Yilmaz term from equation ( |3.2| ) is just given by the 
Newtonian stresses ( |A.5| ) 

t ik = \t$ (4.3) 
so from equation ( A.6| ) one has t jk = XpU^. Also one can evaluate T l ooc 2 = 



Uj from the metric and $ = U/c 2 . The result then is 

p[dv i /dt + v\ k v k } = -P >i -p(l-\)U i . (4.4) 



This is the desired Newtonian result ( |A.3| ) in the Einstein case A = but 
lacks any Newtonian gravitational force in the Yilmaz case A = 1. 



Appendix 

A Newtonian Benchmark 

The Newtonian theory we expect to recover as a limiting case of any rela- 
tivistic theory will be taken to be the hydrodynamics of a perfect fluid. This 
theory has as its gravitational field equation the Poisson equation 

U kk = AttGp (A.l) 

for the generation of gravitational potentials, the continuity equation 

dp/dt + {pv k ) )k = (A.2) 

which shows how fluid motions v k affect the density, and the Euler equation 

thfi/dt = dv l /dt + v\ k v k = -Ui - (l/p)P,i (A.3) 

which shows how the force per unit mass changes the momentum per unit 
mass v l . (Note that for a point mass one has U = —GM/r which is the 
opposite sign from that used for the symbol U in [§, §39.7].) To have a 
complete deterministic theory these equations must be supplemented by an 
equation of state p = p(P). Thermal properties including viscosity and 
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heat transfer are neglected here, and the equation of state is assumed not 
to involve the temperature. 

The Euler equation can be replaced in the Newtonian theory by an equiv- 
alent form in which all forces, not just pressure forces, are described by a 
stress tensor or momentum flux tensor 

2$ = pvV + P» + 4 (A.4) 

where pv % v 3 give the convective transport of momentum and P %3 are the 
mechanical stresses. (When two regions are in contact the forces they exert 
on each other represent, by action equals reaction, a momentum gain for one 
region and a corresponding loss for the other. Thus a contact force per unit 
area is a momentum transfer rate per unit contact area.) The gravitational 
stresses 

4 = (l/ATtG^UiUj - l6 ij U k U k ] (A.5) 

are chosen so that, using the Poisson equation, they give the required grav- 
itational force per unit volume: 

-t*,k = -0-/* irG ) u <i u ,kk = -P u ,i • ( A - 6 ) 

For application to fluids one sets P % i = PS 1 ^ and finds that —P tk ^ = ~~ Pi 
is the pressure gradient force per unit volume. Conservation of momentum 
is then expressed by 

d(pv i )/dt + T$ k = Q (A.7) 

which treats momentum on a "per unit volume" basis in contrast to the 
previous Euler equation flA.3| ) which is on a "per unit mass" basis. It is 
straightforward, using the equation of continuity, to reduce equation ( |A.7| ) 
to equation ( |A.3j ). 
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